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Abstract. We briefly review the Ising model with uncorrelated, quenched random-site or random-bond disorder, which has 
been controversial in both two and four dimensions. In these dimensions, the leading exponent a, which characterizes the 
specific-heat critical behaviour, vanishes and no Harris prediction for the consequences of quenched disorder can be made. 
In the two-dimensional case, the controversy is between the strong universality hypothesis which maintains that the leading 
critical exponents are the same as in the pure case and the weak universality hypothesis, which favours dilution-dependent 
leading critical exponents. Here the random-site version of the model is subject to a finite-size scaling analysis, paying 
special attention to the implications for multiplicative logarithmic corrections. The analysis is fully supportive of the scaling 
relations for logarithmic corrections and of the strong scaling hypothesis in the 2D case. In the four-dimensional case unusual 
corrections to scaling characterize the model, and the precise nature of these corrections has been debated. Progress made in 
determining the correct 4D scenario is outlined. 
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INTRODUCTION 

Continuous phase transitions are characterized by critical exponents. In terms of the reduced temperature t and the 
reduced external magnetic field h, which are measures of how far the system is from its critical point (t,h) — (0,0), 
the standard power-law leading behaviour is as follows: 

specific heat at h = 0: &*>(?) — \t\~ a , (1) 

spontaneous magnetization at h = 0, t < 0: m^t) ~ \t\^ , (2) 
susceptibility at h = 0: %*,{t) — |f|~ y , (3) 

magnetization at t = 0: m^ili) ~ \h\% , (4) 

correlation length at h — 0: £»(0 ~ |*|~ v . (5) 

The subscripts here represent the linear extent of the system. The correlation function at criticality decays as 

^(x)=x- {d - 2+r l\ (6) 

where x represents distance along the lattice, the dimensionality of which is d. These six critical exponents are related 
by four scaling relations, namely 

vd = 2 -a, (7) 
2j5 + Y = 2- a, (8) 
j3(S-l) = y, (9) 
v(2-tj) = y. (10) 



If a pure system (i.e., a system defined on a regular intact lattice) is characterized by a particular set of critical 
exponents, it is interesting to ask what happens to this set if randomization is introduced to the lattice sites or bonds. 



The quenched random removal of sites or the randomization of bond strengths is believed to immitate the presence of 
impurities in real physical systems, so this question is relevant for meaningful comparison with experiments. 

In most cases the Harris criterion provides an answer to this question [1]. If a > in the pure system, quenched 
disorder is relevant and the critical exponents change as such disorder is added. If, on the other hand, a < in 
the pure model, then this type of disorder does not alter critical behaviour and the critical exponents are unchanged 
upon randomizing the lattice structure. The Harris criterion does not, however, provide aa clear an answer in the 
circumstance where specific-heat critical exponent a vanishes in the pure model, as is the case for the Ising model in 
both two and four dimensions. As a result, these circumstances have been quite controversial. In the following two 
sections the histories and natures of the controversies in 2D and 4D are outlined. An approach to tackle these subtle 
issues is presented in the following section where the results of applications of this method to both cases are also 
presented. Conclusions are drawn in the final section. 

THE TWO-DIMENSIONAL CASE 

In the same manner as the pure Ising model in two dimensions can be formulated as a lattice theory of free fermions, 
the randomized model can be formulated as a theory of interacting fermions. This interaction can be considered 
perturbatively about the solvable pure case. In [2] Dotsenko and Dotsenko used a truncated form of Grassmann field 
theory together with the replica trick, to link the problem to the N — Gross-Neveu model. They derived the following 
behaviour for the specific heat in the RBIM: 

c,„(0~ln|ln|f||. (11) 

They also derived an expression for the scaling behaviour of the susceptibility, namely X°°(t) ~ t ~ 2 exp (— c\ In | In \t \ | | 2 ), 
where c is a constant. Their value 7=2 there was different to that in the pure theory (which has 7=7/4) and the 
notion that the leading critical exponents may change when going from the pure to the random 2D Ising model became 
known as the weak universality hypothesis. In fact, the notion of weak universality as advanced by Suzuki [3] is that 
some exponents may change with dilution, but combinations which appear in terms of the correlation length (e.g., in 

X°°{t) ~ %™ V {t )) are dilution independent. I.e., j3 /vand y/v are unchanged (as are <5 and rj). 

In [4], Jug derived two-loop renormalization-group (RG) results in the 2D (and 3D) RSIM. (For a recent review on 
the 3D RSIM, see [5].) He also worked out the exact RG and e-expansion along the curve in (rc,<i)-space where the 
pure system's specific-heat a-exponent vanishes. Here n is the number of components of the order parameter in 0(n) 
models. This includes the n — 1 RSIM case in d — 2 (as well as d — 4, which we shall discuss in the next section). 
For the 2D case, he showed 1 that the critical behaviour is controlled by the pure-model fixed point, so that the leading 
exponents are the same as in the pure case. In particular, and contrary to the susceptibility results of [2], Jug's approach 
gave 7 = 7/4 for the diluted model - which is the same value as in the pure case. The notion that the leading critical 
exponents are unchanged by randomizing the lattice structure became known as the strong universality hypothesis. 

For d — 2, Jug also derived a change in the logarithmic form of the heat capacity, either to | In \t \ \ 01 with a = 0{e) 
or In | In \t || if a = [6]. In [7], Jug confirmed the result (11) for the heat capacity for the 2D RBIM using Grassmann 
field theory. On the basis of the chronology outlined above, we refer to the (now famous) proposed behaviour (1 1) in 
the specific heat as the Dotsenko-Dotsenko-Jug (DDJ) double logarithm. 

Shalaev later introduced bosonization to the above approaches and derived that all critical exponents are the same 
as in the pure case but there are, in fact, non-trivial multiplicative logarithmic corrections to the susceptibility and 
correlation length in the RBIM [8]. He derived the exponents of these logarithms and again obtained the double 
logarithm of the specific heat, results which were later re-obtained by Shankar and Ludwig [9]. Using transfer matrix 
techniques together with the self-duality of the RBIM, and the replica trick, to map the model to the 0(N) Gross-Neveu 
model, to which RG is applied to one loop, Jug and Shalaev later derived the logarithmic corrections for the remaining 
quantities in the RBIM [10]. 

Modifying the standard scaling expressions (l)-(6), above, to include multiplicative logarithmic corrections, we 
write the universal scaling forms 

c(0 ~ |;r a |ln|;|| ft , (12) 



1 In a contribution to the 1983 Geilo School, " Multicritical phenomena" , Jug presented the work of [4] and announced that for the 2D RBIM 
treated with the Grassmann field theory he had reached Co=(f) ~ A In | In \t \ \ +B\ In |f| | _1 with A and B undetermined constants [6]. 



m»(0 ~ I^M*^, (13) 

x ~(t) ~ kriinkir, d4) 

OToo(ft) - h*\ln\h\\ s . (15) 
&.(*) - |f|- v |ln|f|r, (16) 
Sf..(jc) ~ x-( rf - 2+I l)(lnx)V (17) 

Allowing also for the possibility of logarithmic corrections to the FSS behaviour of the correlation length, we also 
write the universal scaling form 

&(0)~L(lnZ,)*. (18) 

A set of universal scaling relations for the correction exponents has recently been developed which connects the 
universal hatted exponents in a manner analogous to (7)-(10). These are [11, 12] 

( 1+dq-dv if a = and ^ n/A 
~ y dq — dv otherwise, 

2j3-y = dq-dv, (20) 

J3(5-1) = 88-Y, (21) 

f) = y-v(2-rj). (22) 

In the first of these, refers to the angle at which the complex-temperature zeros impact onto the real axis. If a = 0, 
and if this impact angle is any value other than n/A, an extra logarithm arises in the specific heat. This is expected to 
happen in d = 2 dimensions, but not in d = 4, where = n/A [12]. In [12] it was also shown that 

if a = and if d(v-q) = l, (23) 

then the specific heat necessarily has the double-logarithmic divergence (11). 

The Jug-Shalaev-Shankar-Ludwig (JSSL) values for the leading critical exponents and for their logarithmic coun- 
terparts are 

a = 0, J3 = i, Y= 7 -, 5 = 15, v = l, n = \, (24) 

« = 0, P = ~^, 7=1' 5=0, v = \, 0=0- (25) 

The observation that the set of leading exponents (24) obeys the usual scaling relations (7)-(10) is a trivial one, since 
they are identical to those of the pure Ising model in d = 2. More interesting is the observation that, with q = 0, 
[13, 14, 15, 16], the JSSL correction exponents (25) obey the scaling relations for logarithmic corrections (19)-(22). 
The observation that they also obey (23) with q = 0, leads to a new route to the derivation of the DDJ double logarithm 
in the specific heat [12]. 

Analytically and numerically based alternatives to the JSSL scenario and to the DDJ double logarithm in the specific 
heat have been made in the literature. Timonin commented that the replica trick, which was employed in the previous 
analytical approaches, encounters a problem in that RG is only valid with 1 or 2 replicas, preventing the necessary 
n — > limit [17]. Using Grassmannian methods and perturbative RG (with certain other assumptions but without 
replicas) he derived a = — 1/2, i.e., a finite specific heat, in the RBIM. Ziegler used a supersymmetric formulation, 
where n = is replaced by N bosons and 2N fermions to give a non-perturbative approach to show that the specific heat 
does not, in fact, diverge in the RSIM [18] and RBIM [19]. Besides the work of Jug, this was the only analytic work 
on the RSIM to this point. Finally, in 1998, Plechko [20] used Grassmann lattice theory to give theoretical support for 
the double logarithm in the RSIM. 

To summarize, both the RSIM and the RBIM have been targeted over many years using a plethora of analytical 
approaches, some of which support the weak universality hypothesis and others of which support strong universality. 
In order to discriminate between the two, and to decide whether or not the specific heat diverges in these models, there 
have also been many numerical investigations of the problem. 

Early numerical work by Zobin was supportive of the strong hypothesis in the RBIM in that the susceptibility 
exponent 7 was found to be unchanged there [21]. In 1990, Andreichenko, Dotsenko, Selke and Wang [22] used such 
a numerical approach to present evidence in support of the strong hypothesis. They focused on the RBIM because 



TABLE 1. Recent works supportive of the weak or strong scaling hypothesis and for and 
against the double logarithm in the specific heat in the RBIM and RSIM in d = 2 dimensions. 



RBIM I RSIM 



Support strong universality hypothesis 


[14, 


16,21,23, 33, 36] 


[37, 38, 


39] 


and theoretical support for a — a = 


[2,"/ 


, 8,9, 10, 12] 


[4, 6, 12, 20] 


or numerical support for a = a = 


[15, 


22, 25, 34, 35] 


[25, 26, 


27] 


Support for weak universality hypothesis 


[40] 




[41] 




and theoretical support for finite C«,(f) 


[19, 


42] 


[18] 




or numerical support for finite G»(f) 


[24, 


28] 


[29, 30, 


31,32] 



the self-duality of that model leads to an exact value for the critical temperature T c , ameliorating some aspects of the 
numerical analyses. Since then many numericists claim support for the strong hypothesis and the double logarithm in 
the specific heat, mostly for the RBIM, but also for the RSIM. However many others support the weak hypothesis and 
a finite specific heat (mostly for the RSIM, but also for the RBIM). The situation is summarized in Table 1. 

While Roder et al. presented strong numerical evidence that f = 7/8 in the RBIM [23], their series-expansion 
approach did not lead to a clear result for the specific heat, and almost any reasonable value a < 1 could be supported 
from their data. Indeed, it was emphasised in [24] that plots of the type contained in Refs. [22, 25] for the RBIM and 
in Refs. [25, 26, 27] for the RSIM, which purport to display double-logarithmic behaviour of the specific heat do not 
actually imply such divergence. I.e., it is very difficult (even impossible) to distinguish between 

Coo(t) ~A + Bln|ln|f|| and c„(t) ~ A -B\t \~ a or c„(0 ~ A -B|ln|?||" , (26) 

with a < or a < 0, on the basis of direct numerical simulations of the specific heat. Numerically based counter-claims 
for the latter two behaviours (so that the specific heat remains finite) in the random-bond [24, 28] and random-site 
models [29, 30, 31, 32] also exist. 

Thus, like the analytical situation, the numerical approaches to the equilibrium lattice-disordered Ising models in 
2D have generated much controversy and debate. (For non-equilibrium scaling aspects in these types of models, see 
[43].) While it is perhaps fair to say that the strong hypothesis is mostly favoured, agreement is not universal. Also, 
the double logarithmic behaviour of the specific heat has resisted attempts at verification because of the difficulties in 
disentangling the scenarios of Eq.(26). 

Here we present a method which circumvents these difficulties. We use the recently developed scaling relations for 
logarithmic corrections [11, 12] to express the scaling of Lee- Yang zeros (in particular their density) in terms of the 
specific heat exponents a and a. Since the density of zeros is not accompanied by a constant or homogeneous term, it 
provides a cleaner Ansatz with which to extract a and a numerically, at least in the d = 2 case. 



THE FOUR-DIMENSIONAL CASE 

Since the upper critical dimensionality of the RSIM, like its pure counterpart, is d = 4, the leading critical exponents 
are given by mean field theory 

a = 0, P = \, 7=1, 8 = 3, v = i, 77=0, A=|, (27) 

and there is no weak universality hypothesis. However, this model is characterised by unusual corrections to scaling 
together with multiplicative logarithmic terms, the precise nature of which are unsettled. 

The consensus in the literature is that the scaling behaviour of the RSIM in four dimensions is given by [4, 44, 45, 
46,48,49] 



CUi) - A — B\t "cx P ; --l x j^-\\n\t\\ )|ln|r|| & , (28) 



*-(0 ~ kr r exp(^|ln| f || |ln| f ||r, (29) 



TABLE 2. Analytic predictions for the logarithmic correction exponents in the 4D RSIM. Entries 
in boldface come directly from the cited references. The remaining entries come from the scaling 
relations for logarithmic corrections. 

Log exponent a /3 y 8 v i) q A 



Aharony [44] 


0.5 


0.25 





0.167 








0.125 


0.25 


Shalaev [8] 


1.237 


0.434 


-0.368 


0.167 


-0.189 


0.009 


0.120 


0.803 


Jug [4] 


0.5 


0.252 


0.005 


0.170 








0.248 


Geldart & De'Bell [46] 


1.246 


0.439 


-0.368 


0.170 


-0.187 


0.005 


0.125 


0.807 


Ballesteros et al [26] 


0.5 


0.255 


0.009 


0.173 





0.009 


0.125 


0.245 



Ut) ~ |fr v exp^A|ln| f ||j|ln| f |r. (30) 
Modification of the theory presented in [1 1] gives that the scaling behaviour for the magnetization in this 4D model is 



m»(0 = f"exp^-ly^Hf||J|lnf|*, (3D 

moo (h) = h^\\nh\ & . (32) 

There is no dispute in the literature regarding the unusual exponential correction terms in (28)-(30), but there are five 
different sets of predictions for the exponents of the logarithmic terms, which differ from their counterparts in the pure 
model. 

Using RG, Aharony derived the unusual exponential terms in (28)-(30), and [44] 

& = i, y=0, v = 0. (33) 

In [45], Shalaev refined Aharony 's calculations to higher order in perturbation theory, yielding 

a = 1.2368, y= -0.3684, f]=0094. (34) 
Jug's calculations along the a = line in (n,d) space yielded [4] 

a =1/2, y = 1/212 w 0.0047 (35) 
at (n,d) = (1,4). In [46], Geldart and De'Bell derived 

a w 1.2463, f« -0.3684, *? = -]- =0.0047, (36) 
and finally Ballesteros et al. [48] gave 

&=\, t=-^~ °- 0094 ' ^ = 0, q= 1 -. (37) 

From these fragmented pictures, complete scaling scenarios may be built using the scaling relations for logarithmic 
corrections. These scenarios are summarized in Table 2. 

It is remarkable that none of the five analyical works on this model agree regarding the detail of the logarithmic 
correction exponents, a circumstance which motivated our investigations into the model. We chose to investigate the 
random-s;fe versions of the models, as this has hitherto proved most controversial in 2D. These investigations use the 
Lee- Yang zeros of the partition function, which we next describe. 



A FRESH PERSPECTIVE: LEE- YANG ZEROS 



We have examined the RSIM in both two and four dimensions from a fresh perspective, namely using the Lee-Yang 
zeros of the partition function. 

A phase transition is a physical manifestation of a mathematical non-analycity in the free energy. Since the free 
energy is essentially the logarithm of the partition function, non-analyticities arise when the latter vanishes. Writing 
the partition function for a system of size L as Zi{t,h), one may consider zeros in either of the variables t or h. In 
each case, the zeros are located in the complex plane. This tack was first suggested by Lee and Yang, and complex-/z 
zeros (which, according to a theorem by the same authors, are usually located along the imaginary /z-axis) are called 
Lee- Yang zeros. Above the critical temperature (f > 0), where there is no phase transition, they are located away from 
the critical point h = 0. There, the zeros may be considered to be "proto-critical points", in the sense that they have the 
potential to become actual transition points. 



Scaling 

The start of the distribution of complex-/; zeros in the t > phase is called the Yang-Lee edge, which we write as 
rYh(t)- As the temperature is reduced towards the critical one (f — > 0), the edge and the distribution of zeros move 
towards the real axis, which they pinch at t = 0. This pinching, which precipitates the phase transition, scales in a 
manner characterised by critical exponents. Allowing for logarithmic corrections here too, we write 

r YL (f)~f A |lnf| A , (38) 

in the 2D case. For the 4D model, we have to account for the unusual correction terms. Following the theory presented 
in [1 1], it turns out that the scaling behaviour for the Yang-Lee edge is (see also [47]) 



r YL (t)~t*exp\--^-\\n\t\\ j in/; 1 (39) 

The exponent A is called the gap exponent. It is related to the other exponents via 

A = J3 + y. (40) 
The logarithmic analogue to this scaling relation is [1 1] 

A = j3-y. (41) 

Besides the scaling behaviour of the Yang-Lee edge, one may also consider the density of zeros, which, for an 
infinitely large system, we write as goo(r, t), where r parametrizes their locus along the imaginary /z-axis (assuming the 
Lee- Yang theorem holds). In fact it is more convenient to consider the integrated, or cumulative, distribution function 
of zeros, which is defined as 

G«,{r,t) = / goo(s,t)ds. (42) 

From [11, 12], in the 2D case at t = this is 

GooM^r^pnr)"- 1 ^ 2 - ^. (43) 
Following a similar approach to that oulined in [1 1], in the 4D case, one determines 



G»(r) ~ r 2 -^ exp M - g U A | i„ r | | i nr |&-(2-«)| (44) 



Using the mean-field values y = 1 and A = 3/2, the exponential term drops out of this expression. 
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FIGURE 1. The leading FSS (left) and corrections (right) for the susceptibility in the RSIM in two dimensions for weak 
(p = 0.88889, upper data set, red online), moderate (p = 0.75, middle set, green online) and strong (p = 0.66661, lower set, 
blue online) dilution values. 



Finite-Size Scaling 

The finite-size scaling (FSS) behavior for the susceptibility, the Yang-Lee edge and the specific heat for each model 
is determined by using (16) and (30) to write the reduced temperature in terms of the correlation length. Then the 
behaviour of the various functions may be expressed directly in terms of %o°(t). For sufficiently large lattices, ^oo(f) 
may be replaced by §l(0) at the critical point. The relationship between §l(0) and the lattice size is, in turn, given by 
(18), where q = in the 2D case and q is listed in Table 2 in the 4D scenarios. 

The unusual exponential correction terms in 4D, which otherwise swamp the logarithmic corrections, are not in 
doubt in the literature. To probe the contested logarithmic correction terms, therefore, these terms have to be removed. 
It turns out that these terms cancel out in the FSS expressions for the susceptibility, the Yang-Lee edge and the density 
of zeros in the 4D case. 

To summarize, the FSS behaviour for the susceptibility and for the edge (i.e., the first Lee- Yang zero) in both models 
is given by 

Jfc(0)~L$|lnL|* 
n {L)~L-$\\nL\P 

The FSS of the specific heat in each model is 

cl(0) ~ (InL)" + constant when d = 2, (47) 
c L (0) « A-B'exp ^-2y^lnLj (InL)" when d = 4. (48) 

In the 2D case, it is difficult to test the Ansatz (47), as explained around Eq.(26). There, we shall instead extract a 
using the form (43) for the density, together with the scaling relations for logarithmic corrections. 



where £ = Yl^Xlll , (45 ) 

„ vA + Av-Aq 
where p = . (46) 



NUMERICAL APPROACH 

We have simulated the RSIM in both d = 2 and d = 4 dimensions using the Wolff single-cluster algorithm [50]. The 
quenched dilution is implemented by occupying a site at a given point on the lattice with probability p, so that p = 1 
corresponds to the pure (intact) model. In the d = 2 case, we simulated at weak, moderate and strong site-dilution 
values, represented by p = 0.88889, p = 0.75 and p = 0.66661, respectively, and we used lattices of size L = 32, 48, 
64, 96, 128, 196 and 256. In the d = 4 case p = 0.8 and p = 0.5 with L = 8, 12, 16, 24, 32 and 48 were used. We also 
simulated the pure p = 1 models. The simulations were carried out at the critical temperatures estimated in [26, 48] 
and periodic boundary conditions were used. Up to 1000 realizations of disorder were generated for each lattice size 
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FIGURE 2. The leading (left) and corrections (right) to FSS for the Yang-Lee edge for the RSIM in 2D for weak (p = 0.88889, 
lower data set, red online), moderate (p = 0.75, middle set, green online) and strong (p = 0.66661, upper set, blue online) dilution 
values. 
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FIGURE 3. The leading (left) and corrections to (right) finite-size scaling for the Yang-Lee edge for the RSIM in two dimensions 
for weak (p = 0.88889, upper data set, red online), moderate (p = 0.75, middle set, green online) and strong (p = 0.66661, lower 
set, blue online) dilution values. 



and each dilution value, with the specific heat, susceptibility, and lowest zeros being determined for each one. These 
were then averaged to give estimates for each quenched model. Further details are supplied in [51, 52]. 



The Two-Dimensional Case 

We begin the analysis with the susceptibility, the FSS for which is given in (45) with the JSSL predictions (24) 
giving y/ v = 7/4. In fact although the weak hypothesis advocates dilution-dependent leading exponents, the ratio y/v 
is believed to be fixed there, so one cannot distinguish between the strong and weak scenarios on this basis. However, 
Roder et al. have presented compelling evidence for the value 7=7/4 [23], which we can input into FSS in order 
to extract v. The results of this process are depicted in Fig. 1 and summarized in Table 3. From fits to the leading 
behaviour Xl{0) ~ lJl v , the theoretical value y/v — 7/4 and hence v = 1 is clearly supported for each dilution value. 
Accepting this value, and fitting for the logarithmic-correction exponent £, one obtains values compatible with zero, 
and therefore compatible with theory. Then inputting the value y = 7/8 clearly evidenced in [23], one obtains from 
(45) and (24), values for v — q clearly compatible with the JSSL value 1 /2 (see Table 3). The scaling relation (19) then 
leads to the estimates a = -0.01 ± 0.04, a = 0.02 ± 0.03, and a = 0.01 ± 0.04, for p = 0.88889, p = 0.75, and for 
p = 0.66661, respectively. 

The FSS for the lowest lying Lee-Yang zeros are plotted in Fig. 2 and the results of the corresponding fits are also 
listed in Table 3. The fits to the leading power-law behaviour are consistent with the theoretical value A/v = 16/8. 
Accepting this leading value, and fitting for the corrections lead to values of p compatible with zero. From Eq.(46), 
then, one can estimate A using the established values of v, A and v — q. These estimates for A are listed in Table 3. 

To summarize the situation at this point, analyses of the FSS of the susceptibility and the Yang-Lee edge yield 



results which are dilution-independent for both the leading and logarithmic-correction exponents and fully compatible 
with the theoretical predictions of JSSL. 

We finally turn to the density of zeros (43). For large enough lattice size, it is expected that the integrated density 
Goo(r) may be estimated from finite lattices as [11] 

G L (rj(L)) = ?tj±, (49) 

where j is the index of the zero r 7 (L). We note that unlike the specific heat, there is no constant term contaminating 
the expression (43), allowing a and a to be cleanly extracted (having already established A and A). 

The density distribution function of zeros is plotted in Fig. 3 using the first four zeros for lattices from size L = 32 to 
L = 256 (28 points in all) for each value of the dilution. Excellent data collapse is evident in each case, and fits indicate 
that each curve goes through the origin, as it should at the critical value of the temperature. The potential logarithmic 
corrections are firstly ignored and fits to the leading scaling of the density are made. For p = 0.88889, fits to the eight 
lowest data points yield and (2 — a)/A = 1 .076(16), corresponding to the estimate a = —0.02(3). The corresponding 
results in the p = 0.75 and p = 0.66661 cases are (2-a)/A= 1.062(10) (a = 0.01(2)) and (2 - a)/ A = 1.066(15) 
(a = 0.00(3)), respectively. These results for a are gathered in Table 3. We next accept the theoretical value a — for 
each dilution and explore potential multiplicative logarithmic corrections by plotting InG — 16/15 In r against In (In r) 
in Fig. 3. The Ansatz (50) becomes 

Goc(r) ~rB(lnr) a . (50) 

A fit to all data points for p = 0.88889 gives a = 0.012(3), four standard deviations from the theoretical value of zero. 
However, focusing on the scaling region closer to the origin establishes compatibility with the DDJ value. For example, 
fitting to the lowest eight data points yields a = —0.02(5). The equivalent results for p = 0.75 and p — 0.66661 are 
a = —0.01(3) and —0.04(5), respectively. The corresponding fits are depicted in Fig. 3 and the estimates for a are 
summarized in Table 3. These values constitute numerical evidence that a = a = 0, independent of dilution and in 
favour of DDJ and JSSL. 



The Four-Dimensional Case 

Starting with the weaker dilution value p — 0.8, ignoring logarithmic corrections, and fitting to the leading form of 
(45) yields the estimate y/v = 2.14 ± 0.01 using lattice sizes L = 8 — 48 (see Fig. 4). Attributing the deviation from the 
mean-field value y/v = 2 to the logarithmic corrections, we find an appropriate fit yields £ = 0.39(3) for 8 < L < 48. 
Thus the FSS logarithmic corrections have moved from the pure value £ = 0.5 towards the theoretical estimates for 
the diluted value, namely £ w 0.25 to 0.26. The same analysis for the FSS of the susceptibility at the stronger dilution 
value p = 0.5 gives similar results: the leading exponent is estimated at y/ v = 2. 1 3 ± 0.02 and the correction exponent 

TABLE 3. Estimates for critical and logarithmic-correction from fits to the scaling behaviour of 
the susceptibility and first Lee- Yang zeros for the 2D RSIM. These estimates agree with the JSSL 
theoretical values. 

Theoretical 



Exponent 


value 


p = 0.88889 


p = 0.75 


p = 0.66661 


r/v 


7/4 = 1.75 


1.747 ±0.007 


1.755 ±0.005 


1.752 ±0.007 


=> V 


1 


1.002 ±0.004 


0.997 ±0.003 


0.999 ±0.004 


I = (vy-yv + rq)/v 





-0.01 ±0.03 


0.02 ±0.03 


0.01 ±0.03 




1/2 


0.51 ±0.02 


0.49 ±0.02 


0.49 ±0.02 


A/v 


15/8 = 1.875 


1.879 ±0.004 


1.878±0.006 


1.878 ±0.006 


p = (vA + Av-A^/v 





-0.01 ±0.02 


0.04 ±0.02 


0.01 ±0.03 


A 


-15/16 = -0.9375 


-0.95 ±0.02 


-0.95 ±0.03 


-0.95 ±0.03 


a 





-0.02±0.03 


0.01 ±0.02 


0.00 ±0.03 


a 





-0.02±0.05 


-0.01 ±0.03 


-0.04 ±0.05 
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FIGURE 4. The leading FSS (left) and corrections (right) for the susceptibility in the RSIM in 4D for weak (p = 0.8, upper data 
set, red online) and strong (p = 0.5, lower set, blue online) dilution values. 



is Q = 0.37(4) for 8 < L < 48. These results are summarised in Table 4, together with results obtained from the same 
fits with the smallest lattices removed. 

The FSS behaviour for the Yang-Lee edge is plotted in Fig. 5. Fitting only to the leading behaviour in (46), for the 
weaker dilution given by p — 0.8, one obtains A/ v = 3 .055 (8 ), using all lattice sizes and at the stronger dilution value 
p = 0.5 we obtained A/v = 3.068(13). Again, we interpret these as supportive of the Gaussian leading behaviour 
y/v = 3 with logarithmic corrections. 

These logarithmic-correction exponents are estimated by fitting to (46), with the various theories in the literature 
indicating that p = —0.125 to —0.13. We find clean evidence in support of this with the estimates p = —0.15(2) and 
p = —0.20(4) at p — 0.8 and p = 0.5 respectively, using L = 8^-8. Dropping the smallest lattice sizes from the FSS 
analysis yields even more convincing results, namely p = —0.17(4) and p = —0.16(5) at weak and strong dilution, 
respectively. Each of these are supportive of [7, 44, 45, 46, 48] and are summarised in Table 4. 

The integrated densities of zeros is calculated in a similar manner to the 2D case and are plotted in Fig. 6 for p = 0.8 
and p = 0.5 alongside the equivalent in the pure p = 1 case. A fit to the leading behaviour G(r) ~ r^ 2_a '/ A yields 
(2 — a)/A= 1.32(3), 1.32(1) and 1.32 (1) for p = 1, p = 0.8, and p — 0.5 respectively. These are compatible with 
the theoretical value 4/3, independent of dilution strength. The errors are too large for us to be confident about the 
equivalent density analysis for the logarithmic corrections and instead we examine the specific heat directly, albeit in 
a rather unusual manner. 

Differentiating Ansatz (48) for the specific heat, one finds its slope vanishes when Cl (0) = A (the asymptote L-*«i) 
and when L = exp ((53/12)a 2 ). The specific heat is plotted for the two dilution values in Fig. 6, from which it is clear 
that the second occurrence of zero slope is for a lattice size smaller than L = 8. On this basis, one may conclude 
a ^ \J 12/53Vln8 s» 0.7, and therefore exclude the values a f=s 1.237 and a s» 1.246 given in [45, 46]. The specific- 
heat curves plotted in Fig. 6 are best fits to the Ansatz (48) with a fixed 1 /2, the alternative value from [4, 44, 48]. 




FIGURE 5. The leading FSS (left) and corrections (right) for the Yang-Lee edge in the RSIM in 4D for weak (p = 0.8, lower 
data set, red online) and strong (p = 0.5, upper set, blue online) dilution values. 
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FIGURE 6. Left: The integrated density of zeros for the pure Ising model (upper curve, red online) and the RSIM in 4D with 
p = 0.8 (middle curve, green online) and p = 0.5 (lower curve, blue online). Right: The specific heat for p = 0.8 (upper curve, red 
onlune) and for p = 0.5 (lower curve, blue online) in the 4D RSIM. 

CONCLUSIONS 

We have presented reviews on the quenched-disordered Ising model in two and four dimensions, where the specific- 
heat exponent of the pure models vanish and no clear Harris prediction for critical behaviour at the phase transitions 
can be made. This circumstance has resulted in both the 2D and 4D models being controversial. 

In the 2D case, the debate has persisted for over thirty years. After confirming the JSSL predictions [4, 8, 9, 10] for 
the logarithmic corrections to scaling for the susceptibility in the random-site version of that model, we determined 
the Lee- Yang zeros to high accuracy and verified that their logarithmic corrections also accord with the JSSL scenario. 

In the 2D model, the precise behavior of the specific heat has been especially controversial and notoriously difficult 
to pin down directly. Using recently developed scaling relations for logarithmic corrections [11, 12], together with 
FSS and Lee-Yang zeros, we have presented an alternative approach, which strongly favours the DDJ scenario 
[2, 4, 6, 7, 8, 9, 10]. These analyses were carried out at weak, moderate and strong dilution values, thereby supporting 
of the strong scaling hypothesis that the exponents are dilution independent. 

In the 4D case, our analysis also strengthens the analytical predictions that the Gaussian fixed point of the pure 
model dominates scaling and that the logarithmic corrections in the RSIM differ from those in the pure model as 
predicted in [4, 44, 45, 46, 48]. Furthermore, we have succeeded in discriminating between some of the detailed 
analytic predictions in the literature, and and our analysis favours the predictions of [4, 44, 48] over those of [45, 46]. 
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TABLE 4. Estimates for FSS exponents for various weak and strong di- 
lution. One expects that Xl ~ L 2 (lnL)^ where £ ~ 0.25 to 0.259 and r\ ~ 
L~^i\nL)P , where p ~ —0.125 to —0.130. (In contrast, the pure p = 1 theory 
is known to have £ = 1/2 and p = — 1 /4.) 



L-range 


8-48 


12-48 


8-48 


12-48 


P 


I 


P 


0.8 


0.39 ±0.03 


0.42 ±0.04 


-0.15 ±0.02 


-0.17±0.04 


0.5 


0.37 ±0.04 


0.40±0.06 


-0.20±0.04 


-0.16±0.05 
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